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Abstract 
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1 Introduction 



The study of Bose-Einstein Condensation (BEC) in random media has been an important 
area for a long time, starting with the papers by Kac and Luttinger, see [T] , [2] , and then by 
Luttinger and Sy [3]. In the last reference, the authors studied a non-interacting (perfect) one 
dimensional system with point impurities distributed according to the Poisson law, the so- 
called Luttinger-Sy model. The authors conjectured a macroscopic occupation of the random 
ground state, but this was not rigorously proved until [S]. Although the free Bose gas (i.e., 
the perfect gas without external potential) does not exhibit BEC for dimension less than 
three, the randomness can enhance BEC even in one dimension, see e.g. [I]. This striking 
phenomenon is a consequence of the exponential decay of the one particle density of states 
at the bottom of the spectrum, known as Lifshitz tail , or "doublelogarithmic" asymptotics, 
which is generally believed to be associated with the existence of localized eigenstates [TT] . 

BEC, however, is usually associated with a macroscopic occupation of the lowest one- 
particle kinetic-energy eigenstates, which are spatially extended (plane waves). Therefore, 
it is not immediately clear whether the phenomenon discovered in random boson gases, i.e. 
macroscopic occupations of localized one-particle states, has any relation to the standard 



2 



BEC. This is of particular interest in view of the well-known Bogoliubov approximation 
[H] and its applications to disordered boson systems, see e.g. [7J, [8], where the a priori 
assumption of the momentum-space condensation is essential, and is far from trivial to 
check. 

In this paper, we prove that for the perfect Bose gas in a general class of non-negative 
random potentials, BEC in the random localized one-particle states and BEC in the lowest 
one-particle kinetic-energy states occur simultaneously, and moreover the density of the 
condensate fractions are equal. Our line of reasoning is also applicable to some non-random 
systems, for example to the case of the perfect gas in weak (scaled) external potentials 
studied in [T7j . 

The structure of the paper is as follows: in Section [5] we describe our disordered system, 
and in Section [31 we recall standard results about the corresponding perfect Bose gas. The 
existence of generalized BEC in the eigenstates of the one-particle Schrodinger operator 
follows from the finite value of the critical density for any dimension, which is a consequence 
of the Lifshitz tail in the limiting Integrated Density of States (IDS). It is well-known that 
the IDS is a non-random quantity, see e.g. [TT], and therefore the BEC density is also non- 
random in the thermodynamic limit. In Section HJ we turn to the main result of this paper: 
we show that this phenomenon occurs if and only if there is also occupation of the lowest one- 
particle kinetic-energy eigenstates. The latter corresponds to the usual generalized BEC in 
the free Bose gas, that is a perfect gas without external potential. To establish this we prove 
the existence of a non-random limiting occupation measure for kinetic energy eigenstates, 
and moreover, we obtain an explicit expression for it. To this end, we need some estimates for 
the IDS before the thermodynamic limit, namely a finite volume version of the Lifshitz tail 
estimates, which we prove in Section [5J using techniques developed in [9], \IU\. For any finite 
but large enough system, these bounds hold almost surely with respect to random potential 
realizations. In Section El we look at the particular case of the Luttinger-Sy model and 
examine the nature of the condensate in the one-particle kinetic energy eigenstates, showing 
that although there is generalized BEC, no condensation occurs in any of them. In Section [7J 
we describe briefly how the method developed in Section H] applies with minor modifications 
to a perfect Bose gas in a general class of weak (scaled), non-random external potentials. 
To make the paper more accessible and easy to read, we postpone some technical estimates 
concerning random potentials and Brownian motion to Appendices A and B, respectively. 



2 Model, notations and definitions 

Let {A; := (— 1/2, l/2) d }i^ 1 be a sequence of hypercubes of side I in M. d , d ^ 1, centered at 
the origin of coordinates with volumes Vi = l d . We consider a system of identical bosons, 
of mass m, contained in A;. For simplicity, we use a system of units such that h = m = 1. 
First we define the self-adjoint one-particle kinetic-energy operator of our system by: 

h° := -IA D , (2.1) 

acting in the Hilbert space M{ := L 2 (Ai). The subscript D stands for Dirichlet boundary 
conditions. We denote by {Vi;4}fc^i the set of normalized eigenf unctions and eigenvalues 
corresponding to h®. By convention, we order the eigenvalues (counting the multiplicity) as 

£ U 4 ^ 4 • • • • 
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We define an external random potential : f2 x W 1 — > R, x i— > v w (x) as a random 

field on a probability space (fi, T, P), satisfying the following conditions: 

(i) v w ,u G f2, is non-negative; 

(ii) p := F{uj : i/"(0) = 0} < 1. 

As usual, we assume that this field is regular, homogeneous and ergodic. These technical con- 
ditions are made more explicit in Appendix [B] Then the corresponding random Schrodinger 
operator acting in J$? := L 2 (M. d ) is a perturbation of the kinetic-energy operator: 

h u := -|A + v", (2.2) 

defined as a sum in the quadratic-forms sense. The restriction to the box A;, is specified 
by the Dirichlet boundary conditions and for regular potentials one gets the self-adjoint 
operator: 

h? := (-i A + v») D = h? + v\ (2.3) 

acting in M[. We denote by {4%'\ Ef' 1 }^ the set of normalized eigenfunctions and corre- 
sponding eigenvalues of hi. Again, we order the eigenvalues (counting the multiplicity) so 
that Ei < E%' E% h .... Note that the non-negativity of the random potential implies 
that E%' > 0. So, for convenience we assume also that in the thermodynamic limit almost 
surely (a.s.) with respect to the probability P, the lowest edge of this random one-particle 
spectrum is: 

(iii) a.s.-lim^oo E^ = 0. 

When no confusion arises, we shall omit the explicit mention of / and uj dependence. Note 
that the non-negativity of the potential implies that: 

(a) Q{hf) C Q(h°), Q being the quadratic form domain, (2.4) 

(b) fahfip) > (<p,tit<p),V<peQ(hf). 

Now, we turn to the many-body problem. Let := j^j(^f) be the symmetric Fock space 
constructed over M{. Then Hi := dT(hf) denotes the second quantization of the one-particle 
Schrodinger operator hf in &\. Note that the operator Hi acting in has the form: 

Hi = Y,K l (2-5) 

where a*(0j),a(0j) are the creation and annihilation operators (satisfying the boson Canon- 
ical Commutation Relations) in the one-particle eigenstates {4>i := 0^''}i>i of hf. Then, the 
grand-canonical Hamiltonian of the perfect Bose gas in a random external potential is given 
by: 

Htiji) := Hi - uNi = Y,{ E f - Wi) (2-6) 

where iVj(<^) := a*(0j)a(0,) is the operator for the number of particles in the eigenstate (pi, 
Ni := Yli Ni((f>i) is the operator for the total number of particles in A; and \x is the chemical 
potential. Note that Ni can be expanded over any basis in the space M{, and in particular 
over the one defined by the free one-particle kinetic-energy eigenstates {tp l k ,Sk}k- 



1 



We denote by (— )/ the equilibrium state defined by the Hamiltonian Hi(p) 



Tr^ exp(-/3^(^)) 

For simplicity, we shall omit in the following the explicit mention of the dependence on the 
thermodynamic parameters (/3,p). Finally, we define the Thermodynamic Limit (TL) as the 
limit, when I — > oo. 



3 Generalized BEC in one-particle random eigenstates 

In this section we consider the possibility of macroscopic occupation of the one-particle ran- 
dom Schrodinger operator (12.31) eigenstates {</>«}«>!. Recall that the corresponding limiting 
IDS, v(E), is defined as: 

u(E) := lim vf(E) = lim U{i : Ef <: E}. (3.1) 

t— >oo / — >oo Vi 

Although the finite- volume IDS, if(E), are random measures, one can check that for ho- 
mogeneous ergodic random potentials the limit (13. ip has the property of self- averaging [TTj . 
This means that v(E) is almost surely (a.s.) a non-random measure. Let us define a (ran- 
dom) particle density occupation measures mi by: 



mi (A) := 1 W*))*, ACR. (3.2) 



Vi 

i:E t <=A 



Then using standard methods, one can prove that this sequence of measures has (a.s.) a 
non-random weak-limit m, see (13.81) below. Moreover, if the critical density 

is finite, then one obtains a generalized Bose-Einstein condensation (g-BEC) in the sense 
that this measure m has an atom at the bottom of the spectrum of the random Schrodinger 
operator, which by (iii), Section [21 is assumed to be at 0: 

m({0}) = lim lim £ ± = { ° j2 > ^ (3-4) 

5J.0 «^oo . ^— ^ v ; [ p — p c n p ^ p c 

where p denotes a (fixed) mean density [I], [5]. Physically, this corresponds to the macro- 
scopic occupation of the set of eigenstates 0j with energy close to the ground state 0i. 
However, we have to stress that BEC in this sense does not necessarily imply a macroscopic 
occupation of the ground state. In fact, the condensate can be spread over many (and even 
infinitely many) states. 

These various situations correspond to classification of the g-BEC on the types I, II and III, 
introduced in eighties by van den Berg- Lewis- Pule, see e.g. [12] or [U], [13] • The most striking 
case is type III when generalized BEC occurs in the sense of equation (13 .4p even though none 
of the eigenstates 4>i are macroscopically occupied. The realization of different types depends 



5 



on how the relative gaps between the eigenvalues Ei at the bottom of the spectrum vanishes 
in the TL. To our knowledge, analysis of this behaviour in random system has only been 
realised in some particular cases, see [3J for a comprehensive presentation. The concept of 
generalized BEC is more stable then the standard one-mode BEC, since it depends on the 
global low-energy behaviour of the density of states, especially on its ability to make the 
critical density ( 13. 3p finite. We note also that, since the IDS ( 13.11) is not random, the same 
it true for the amount of the g-BEC (13 ,4p . 

We can also obtain an explicit expression for the limiting measure m. Note that we have 
fixed the mean density p, which implies that we require the chemical potential p to satisfy 
the equation: 

- p = {N , )liM = l^—^l— , (3.5) 

for any I. Since the system is disordered, the unique solution pf := pf(j3,p) of this equation 
is a random variable, which is a.s. non-random in the TL [I], [3J. In the rest of this paper 
we denote the non-random p^ := a.s.-linij_ >0O /Xj'. By condition (iii), Section[21 and by (13.71) 
it is a continuous function of p : 

MM = i- n ^-- PC ' ( 3 - 6 ) 
H-oo\H,F) \ ^ < if p < p c , K J 

where p := p(/3, p) is a (unique) solution of the equation: 

~ P = r e**-M) - 1 ^ ' (3 - 7) 

for p < p c . 



Remark 3.1 Note that p^ is non-positive $3. 6]) . which is not true in general for the random 
finite-volume solution pf. Indeed, the only restriction we have is that pf < E^ 1 , which is 
the well-known condition for the pressure of the perfect Bose gas to exist. We return to 
this question in Section^for the case of random BEC in the free one-particle kinetic- energy 
operator eig en functions. 



We also recall that for (13.61) the explicit expression of the weak limit for the general particle 
density occupation measure is: 

mm = ( «v^? (d f ? t { ;Zr 1,-1 u{AE) ?- >p " (3.s) 

v ' 1 {eP^ i* 00 ' - ly 1 u(dE) if p < p c . v ' 



We end this section with a comment about the difference between the model of the perfect 
Bose gas embedded into a random potential and the free Bose gas. In the latter case, one 
should consider the IDS of the one-particle kinetic-energy operator ( 12.11) . which is given by 
the Weyl formula: 

v°{E) = C d E d l\ (3.9) 

where is Cd is a constant term depending only on the dimensionality d. It is known that for 
this IDS, the critical density (13. 3p is finite only when d > 2, and hence the fact that BEC 
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does not occur for low dimensions. On the other hand, a common feature of Schrodinger 
operators with regular, stationary, non-negative ergodic random potentials is the so-called 
Lifshitz tails behaviour of the IDS near the bottom of the spectrum. When the lower edge 
of the spectrum coincides with E = (condition (hi)), this means roughly that (see for 
example [TTj): 

u(E) ~ e - a/Ed/2 (3.10) 

for small E and a > 0. Hence, the critical density (13.31) is hnite in any dimension, and 
therefore enhances BEC in the sense of (13 .411 even for d — 1,2. This was shown in [I], [5], 
where some specific examples of one- dimensional Poisson disordered systems exhibiting g- 
BEC in the sense of (13.41) were studied. In this article we require only the following rigorous 
upper estimate: 

lim (-E d/2 ) hx(u(E)) > a > , (3.11) 

E-+0+ 

for some constant a. This can be proved (see [9]) under the technical conditions detailed in 
Appendix [Bj which are assumed throughout this paper. In particular these conditions are 
satisfied in the case of Poisson random potentials with sufficiently fast decay of the potential 
around each impurity. 



4 Generalized BEC in one-particle kinetic energy eigen- 
states 

4.1 Occupation measure for one-particle kinetic energy eigenstates 

Similar to (13. 2j) . we introduce the sequence of particle occupation measure rhi for kinetic 
energy eigenf unctions {ipk '■= ip[}k£A*'- 

™^ A ) := V ( N iM)i , Act, (4.1) 

but now in the random equilibrium states (— )/ corresponding to the perfect boson gas with 
Hamiltonian (12. 5p . 

Note that, contrary to the last section, the standard arguments used to prove the existence 
of a limiting measure in TL are not valid for (14.11) . since the kinetic energy operator (12.11) 
and the random Schrodinger operator (12. 3ft do not commute. 

We remark also that even if we know that the measure m (13.81) has an atom at the edge 
of the spectrum (g-BEC), we cannot deduce that the limiting measure m (assuming that it 
exists) also manifests g-BEC in the free kinetic energy eigenstates 

Our motivation to study this problem is that in view of the well-known Bogoliubov ap- 
proximation, it is the later that is required, see [6]- [8]. Indeed, let the second- quantized form 
of the interaction term is expressed in terms of creation/annihilation operators in states ip l k , 
eigenfunctions of kinetic-energy operator (12. ip . Then the so-called first Bogoliubov approxi- 
mation (Bogoliubov ansatz) assumes that only terms involving creation and/or annihilation 
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operators of particles in the ground state ip[ are relevant. The Bogoliubov theory is nontriv- 
ial if there is macroscopic occupation of this zero-mode kinetic-energy operator ground-state. 
Therefore, the g-BEC in the sense (13 .4p is not sufficient to apply the Bogoliubov ansatz. 

Now we formulate the main result of this section. Let 

fif^ == {£ = e(0) = x, £(T) = x'} 

be the set of continuous trajectories (paths) {£(s)}J =0 in M. d , connecting the points x, x', 
and let w T denote the normalized Wiener measure on this set. 



Theorem 4.1 The sequence of measures mi converges a.s. in a weak sense to a non-random 
measure rh, which is given by: 

= f (p - pMde) + F(e)de ifp > Pc 
v ' \ F(e)de if p < p c 

with density F(e) defined by: 

F{e) = (26)^ [ dag(V2ln a ) . 
Jsi 



Here, S\ denotes the unit sphere in M, d centered at the origin, n a the unit outward drawn 
normal vector, and da the surface measure of . The function g is as follows 

1 f 1 1 a; II 2 (1 /2n/3) r 

**> = L Axe ' k * S^^M / n;ti ^(d0e-//-<«.»},(4.2) 

with expectation E> w on the probability space (fi,jF, P). 



Note that since the measures w n ^ on ^q x \ are normalized, we recover from ( 14. 2 p the 
expression for the free Bose gas if we put t> w = 0. 

Before proceeding with the proof, we give some comments about these results. 

(a) First, the existence of a non-trivial limiting kinetic energy states occupation measure 
provides a rigorous basis for discussing the macroscopic occupation of the free Bose gas 
eigenstates. 

(b) Moreover, both occupation measures ( I3.8P and ( 14. ip do not only exhibit simultaneously 
an atom at the bottom of the spectrum, but these atoms have the same non-random weights. 
It is quite surprising that the generalized BEC triggered by the Lifshitz tail in a low dimension 
disordered system produces the same value of the generalized BEC in the lowest free kinetic 
energy states. 



4.2 Proofs 



We start by expanding the measure m in terms of the random equilibrium mean-values of 
occupation numbers in the corresponding eigenstates 0j. Using the linearity (respectively 
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conjugate linearity) of the creation and annihilation operators one obtains: 



rhi(A) = — 



(4.3) 



= v Yl (fa'&Mi, M (a*(<f>i)a((f>j))i 



i,j k:e k eA 



= ^EE \(Mk)\ 2 (a*(<j>Mfa))i. 



In the last equality, we have used the fact that [Hi(fj) , Ni(fa)] = for all i, which implies 
that: 



This is the analogue of the momentum conservation law in the free Bose gas. Although, it 
has a different physical meaning: the conservation of the particle number in each random 
eigenstate fa. 

We first prove two important lemmas. In neither of them we shall assume that the sequence 
fhi has a weak limit, instead we consider only some convergent subsequence. Note that at 
least one such subsequence always exists, see [T6] . 

The first result states that if there is condensation in the lowest random eigenstates {4>i}i, 
then there is also condensation in the lowest kinetic-energy states {ipk}k- Moreover, the 
amount of the latter condensate density has to be not less than the former. 

Lemma 4.1 Let {rhi r } r >i be a convergent subsequence. We denote by m its (weak) limit. 



Proof: Let 7 > 0. Using the expansion of the functions ipk in the basis {(pi}i>i , we obtain: 



(a*(fa)a(^j))i 



tii^j. 



Then: 








for any 5 > 0. The non-negativity of the random potential (12.41) implies: 




|(0,,^)| 2 ^ V ^|(^,^)| 2 ^IVe fe |(^,^)r = -(^,tt)^ -(4>i,hffa) 



1 
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We then obtain: 



™([0,7]) > lim ±- E (^(^(1- E K^*)D 

^ Um E Wr(&)>lr [t-Ei/l) 

^ lim (1- 5h) i- E = (l-*/7MM)>0. 

But 5 is arbitrary, and the lemma follows by letting 5^0. □ 

In the next lemma, we show that the measure rh (14. II) can have an atom only at zero 
kinetic energy 

Lemma 4.2 Let {fhi r } r >i be a convergent subsequence, and rh be its (weak) limit. Then, it 
is absolutely continuous on M + := (0, oo). 

Proof: Let A to be a Borel subset of (0, oo), with Lebesgue measure 0, and let a be such 
that inf A > a > 0. Then: 

m r (A) = TT E ( N lrM)lr (4-4) 
lr k:e k &A 



= t E Ei(^)i 2 w<^ 

lr k:e k &A i 

= v~ S E i(^^)i 2 (iva^)k + ^- e E i(0^fc)i 2 (^,(</>*)k. 

' r fc:e fc eA i:Ei?^a lr k:e k &A i:Ei>a 

for some a > 0. Next, we use (12. 4ft to get the following estimate: 



Ei = {<f>i,hf<f>i) > {fahfyi) = E £ fcl(^^)| 2 ^ a E K<^ 



)\ 2 



Since the the equilibrium value of the occupation numbers {Ni(4>i))i = {e Ei ^ — 1} 1 are 
decreasing with i, the estimate (14.41) implies: 



< ±r- E EiWUMir + (NUK))ir^r E 1 > ( 4 _5) 

lr j:£ 4 ^Q lr k:e k &A 

where <f)i a denotes the eigenstate of hf with the smallest eigenvalue greater than a. Using 
again the monotonicity and the finite-volume IDS (13. lj) we can get an upper bound for the 
mean occupation number in the second term of (14.51) . since: 

p = > y E w&)>' > WiJh ^» • (4.6) 

i i:Ei^.a 

Combining (14.51) and (14. 6 j) we obtain: 



« ^;(«) j a 
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Since the measure v (13.91) is absolutely continuous with respect to the Lebesgue measure, 
and v{a) is strictly positive for any a > the limit r —> oo in (14.71) gives: 

m(A) < — , 
a 

But a > can be chosen arbitrary small and thus rh(A) = 0. To finish the proof, note that 
any Borel subset of (0, oo) can be expressed as a countable union of disjoint subsets with 
non-zero infimum. Our arguments than can be applied to each of them. □ 



Remark 4.1 Lemmas \4-l\ and \4-2\ are fairly general, since they require only the non-negativity 
of the potential and in the random case, ergodicity. In particular they apply to non-random 
(weak) external potentials, that we consider below, as well as to some models of an inter- 
acting Bose gas, as long as the many-particles Hamiltonian still satisfies the commutation 
relation [Hi(fi), Ni((j>i)] = 0. In particular this holds in the case of models with interactions, 
which are diagonal in the occupation number operators. 



Above we exploited the fact that the sequence {rhi}i>i has at least one accumulation point. 
However, to prove convergence, we need to make use of some particular and explicit features 
of the perfect Bose gas, as well as more detailed information about the properties of the 
external (random) potential. In particular, we shall need some estimates of the (random) 
finite volume integrated density of states, see Lemma 15.11 

To this end let us denote by Pa the orthogonal projection onto the subspace spanned by 
the one-particle kinetic energy states ipk with kinetic energy e(k) in the set A. Then using 
the explicit expression for the mean occupation (a*((p i )a(4>i))i and (14. 3 j) we obtain: 

mi{A) = -y Tr Pa ( e W"«) - l)" 1 = V A Tr P A ( e -»W-«)) . (4.8) 

Now we split the measure (14.81) into two parts: 

fh\ = fhi + fh!f^ with , (4.9) 
mf\A) := ViTrP.fe-'W-^l^a/n), 

mf\A) := Y A Tr P A ( e -*W-«)) 1(^ > l/ n ) . 

Note that since the chemical potential satisfies equation (13.51) . fii := fj,f, therefore the indi- 
cator functions l(/i/ ^ l/n) and l(/i/ > l/n) split the range of n into the sums (14. 9ft in a 
random and volume-dependent way. 

We start with the proof of existence of a weak limit of the sequence of random measures 

~ (i) 



Theorem 4.2 Let random potential satisfy the assumptions (i)-(iii) of Section 2. Then 
for any d ^ 1, the sequence of Laplace transforms of the measures fh\ : 

fifcPtm) := / m, (1) (d5) e- te (4.10) 
Jr 
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converges for any t > to a (non-random) limit f(t; /?, /^oo) , which is given by: 

(4.11) 

i/ere denotes the expectation with respect to realizations ( configurations ) u of the random 
potential. Note that the sum on the right-hand side converges for all (non-random) ^ < 0, 
including 0, which corresponds to the case p ^ p c . 



Proof: By definition of Pa the Laplace transformation (I4.10p can be written as: 

/,(*; P,n) = V ±Tre*t ( e -»W-«>) 1(^ <C l/n) . (4.12) 
% Vl 

Now we have to show the uniform convergence of the sum over n to be able to take the term 
by term limit with respect to /. Since for any bounded operator A and for any trace-class 
non- negative operator B one has Tr AB ^ ||v4|| TiB, we get 

ai(n) := ^Tr e~ th ° e^^-^Kjn < l/n) (4.13) 
<: l-Tre-^W-rtlfa ^ l/n) . 

For ~p < p c , the uniform convergence in (14.111) is immediate. Indeed, for I large enough, 
the chemical potential satisfies \i\ < /ioo/2 < 0, which by (13.11) provides the following a.s. 
estimate for (14. 13ft : 

ai{n) < e n/W2 I vf(dE) e~ pE < K x e n ^ /2 , (4.14) 

i[0,oo) 

with some constant K\. 

However, for the case p ^ p c , this approach does not work, since, in fact, for any finite / 
the solutions \ii = pjf of equation (I3.5P could be positive with some probability, event though 
by condition (iii) (see Section 2) it has to vanish a.s. in the TL. We use, therefore, the 
bound: 

ai(n) ^ a i( n ) + a K n ) 5 

a}{n) := 1^ £ e~^\ 

a?(n) := \f £ e~^\ 

which follows, for some < rj < 1, from the constraint ^ 1 due to the indicator function 
in ( 14.13p . Then the first term is bounded from above by: 

a}(n) ^ e?vf (n^ 1 ) . 
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On the other hand, by Theorem 15.11 (finite-volume Lifshitz tails), for a > and < 7 < 
d/2, there exists a subset Q C f2 of full measure, P(fi) = 1, such that for any uj G Q there 
exists a positive finite energy £{u) := £ arj {uj) > for which one obtains: 



for all < £ . Therefore, for any configuration w G (i.e. almost surely) we have the 
volume independent estimate for all n > Af(u) := Eioj) 1 '^ 1 *: 

a]{n) ^ e^e- an(1 "" )7 . (4.15) 

To estimate the coefficients a 2 {n) from above , we use the upper bound: 

a 2 (n) ^ [ v?(dE) e~ nf3E ^ e^' 2 [ v?(dE) e~ n/3E / 2 

J[l/n 1 ~v t oo) J [l/n 1 -^ ,00) 

J [0,00) 

Then for some K 2 > independent of / we obtain: 

a 2 (n) < Kte-W 2 . (4.16) 

Therefore, by (I4.14p in the case p < p c , and by f)4. 15j) . (I4.16P for p > p c , we find that there 
exists a sequence a(n) (independent of I) such that: 

a/(n) ^ a(n) and ^^a(n) < 00 . (4.17) 

Thus, the series (I4.12p is uniformly convergent in I, and one can exchange sum and the limit: 

00 00 
lim fi(t) = lim \ a>i{n) = / lini ai(n) . 



I— >oo I— too ' * ' J I— too 

n=0 n=0 



The rest of the proof is largely inspired by the paper [I]. Using the Feynman-Kac repre- 
sentation, we obtain the following limit: 

lima,(n) = lim -^Tr e~ th ° e - n/3( ^~ w) l(/ii ^ l/n) (4.18) 

l— too l^oo V\ 

= lim \- f [ dxdx'e- th °(x,x')e- n ^ h i-^(x',x) 
1^00 Vi J Al J Al 



1 

'112/ 



Iff -\\x-x'\\ 2 (l/2nf3+l/2t) 

e n ^°° lim — / / dxdx' 



1^00 V, J Ai J A[ (4nHnf3y/2 
x / ^(dOe-^ d -^«W) XAiin/3 (0 f w'idftxAtA?) 



where we denote by Xa ( ,t(0 the characteristic function of paths £ such that G A; for 
all < t < T. Using Lemma IA.2[ we can eliminate these restrictions, and also extend one 
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spatial integration over the whole space: 

lim ^-Tre- th ° e ""W-w) = (4.19) 



V, 



r 1 C p ~\\x-x'\\*(L/2np+L/2t) r 

e^-lim dx- dx' , - — m / w n Hdi)e-^ "^^KW) . 

(x f ,x) 

Now, by the ergodic theorem, we obtain: 

lim ai{n) = lim -Tre~ th ° e - nf> W-n) (4.20) 

l— >oo l-^oo vi 

= e^lim-"- / dx'M dx 6 - - — / ^(de)e-^d- w «W) 

^°°VJ A , }J Rd (Stripy/* J n n, 



e 



x|| 2 (l/2n/3+l/2t) 



(4^)^ L f w We / 



We then get the explicit expression for the limiting Laplace transform: 

which finishes the proof. □ 

Corollary 4.1 For any~p the sequence of random measures rhf^ converges a.s. in the weak 
sense to a bounded, absolutely continuous non-random measure rh^ 1 ', with density F(E) given 
by 



F(e) = (2e) d / 2 - 1 / dag(V*n a ). 
Jsl 



Here, S} denotes the unit sphere in M. d , n a the outward drawn normal unit vector, da the 
surface measure on and the function g has the form 



I C ~ 1 1 *^ 1 1 ^ / ^^J^ /* 

= kw I Axe ' k ' g"* kw* 4 ■""" (df)e " r <"" {m Y^ 



Proof: By Theorem 14.21 the existence of the weak limit fh> l > follows from the existence of 
the limiting Laplace transform. Moreover, we have the following explicit expression: 

f f — \\x\\^ /2t — \\x\\^ /2n(3 I* 

/ rh^(de) e" te = / dx % ~ Ve* \ -— - eJ / w n ^)e~^ d -'^ s ») 



dre -t\\r\\ 2 /2 r d-l I dag ( rn ^ 
[0,oo) JS\ 

= [ de e~ t£ (2e) d/2 ~ 1 [ da g{V2~e n a ) , 

</[0,oo) JSl 

which proves the corollary . □ 
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Corollary 4.2 The measure rrS 1 ^ satisfies the following property: 

7[0,oo) I Pc tfp>Pc 

Proof: By virtue of (14.121) we have: 

/ m«(de) = m^tioo) = Hm V 1 -Re-nW-w) 1(//J ^ i/ n ) . 



Note that by uniformity of convergence of the sum, see (14.151) . (I4.16p . we can take the limit 
term by term (for any value of p), and then: 



/ m (1) (de) = V lim - Tr e -»W-w) = 

[ u(dE) e-^-^ = f u(dE) (e^ E ~^ - I)' 1 

</[0,oo) J\0,oo) 



n2X •'[0.°°) J[0,oo) 

where we use Fubini's theorem for the last step. □ 

We are now ready for the proof of the main result of this section: 

Proof of Theorem I4.ll We first treat the case p < p c . In this situation, the measure mf^ 
is equal to for / large enough, see (I4.9I) . since the solution Mmi^^fif of the equation (I3.5I) 
in the TL is a.s. strictly negative. Thus, the 1 
and the Theorem follows from Corollary 14. 1[ 



in the TL is a.s. strictly negative. Thus, the total occupation measure mi is reduced to m| 



Now, consider the case p ^ p c . Choose a subsequence l r such that the total kinetic-energy 
states occupation measures rhi r converge weakly and a.s., and let the measure m be its limit. 
By Corollary 14.11 all subsequences of measures m, converge to the limiting measure fhP^ . 



Therefore, by 1 14. 9ft . we obtain the weak a.s. convergence: 



lim =: 



By Lemma l4.2[ we know that the measure m is absolutely continuous on (0, oo), and by 
Corollary 14.11 that fhS l > is absolutely continuous on [0, oo). Therefore we get: 

m a - c - = ttjW + m (2)ax - , 

where a.c. denotes the absolute continuous components. 

By definition of the total measure (I4.9p . m([0, oo)) = p and by Lemma |4~T| m({0}) > p— p c . 
Thus, m((0, oo)) < p c and by Corollary 14.21 we can then deduce that the measure has 
no absolutely continuous component and therefore consists at most of an atom at e = 0. 
Consequently, the full measure m can be expressed as: 

m = m a,c ' + bS = fhS^ + bS , 

and since by Corollary 14.21 

b = p - ml c -(de) = p - / m\l\de) = p - p c 
Jm + Jr + 
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for the converging subsequence rhi rJ we have: 

lim rhi r = m (1) + (p - p c )5 . 

By (14.221) and Corollary 14.11 this limit is independent of the subsequence. Then, the limit 
of any convergent subsequence is the same, and therefore, using Feller's selection theorem, 
see [H], the total kinetic states occupation measures rhi converge weakly to this limit. □ 



5 Finite volume Lifshitz tails 

In this section, we give the proof of one important building block of our analysis, Theorem 
15.11 about the finite-volume Lifshitz tails. Recall that this behaviour is a well-known feature 
of disordered systems, essentially meaning that for Shrodinger operators which are semi- 
bounded from below, there are exponentially few eigenstates with energy close to the bottom 
of the spectrum. To our knowledge, however, this is always shown only in the infinite- 
volume limit, see e.g. [TT]. Here, we derive a finite-volume estimate for the density of states, 
uniformly in /, though it could be trivial for small volumes. As one would expect our result 
is weaker than the asymptotic one, in the sense that we prove it for Lifshitz exponent smaller 
than the limiting one. 

Theorem 5.1 Let the random potential v w satisfy the assumptions (i)-(iii) of Section 2. 
Then for any a > and < 7 < d/2, there exists a set Q C Q of full measure, ¥(Q) = 1, 
such that for any configuration uj G Q one can find a positive finite energy S(u) := S an (uS), 
for which one has the estimate: 

v?{E) ^ e~ a l E1 

for all E < S(uj) and for all I. 

For the proof, we first need a result from [9]. 

Lemma 5.1 By assumption (ii) [Section 2) one has, 

p = Pjcu : u w (0) = 0} < 1. 

Let a > p/(l — p), B — 7r/(l + a), and Ei' 1, := E^' be the first eigenvalue of the random 
Schrddinger operator 112. ty) with Neumann (instead of Dirichlet ) boundary conditions. Then, 
for I large enough, there exists an independent of I constant A = A(a), such that 

f{uj:E^ n <B/1 2 } < e~ AV ' . (5.1) 

Detailed conditions on the random potential and a sketch of the proof of this lemma are 
given in Appendix [B] Now we use Lemma I57T1 to prove the following result: 
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Lemma 5.2 Assume that the random potential satisfies the assumptions of Lemma \5.1\ 
Then for any a > and < 7 < d/2, 

^P<{jjz : Ef < 1/nJ > Vie~ an '' , for some I ^ l\ < 00 . 



Proof: Notice that 



F : ^ < I/"} > Vi e~ anl , for some / > l} = ^ P I \J S? \ , (5.2) 



where S 1 / 1 is the set 



i:Ef <^\>V l e- 



The sum in the right-hand side of (15.21) does not provide a very useful upper bound, since 
the sets S] 1 are highly overlapping. We thus need to define a new refined family of sets to 
avoid this difficulty 

To this end we let [a} + be the smallest integer ^ a, and we define the family of sets: 

{( Lo,\{ke anl ) 1 / d \ 1 1 

u: tt{i:^ 1 < nJ >H 

Let fc := [V/e _cm7 ] + . Since Vi = l d , this implies that hf ^ ^ fceCm 7)i/ d j > an d therefore: 

z:£- 4 J+ <-l > Mi:E?> l <- 



n J 1 n 
If now uj G 5", then by the definition of k we obtain: 



% : Ef' 1 < - \ 



since the left-hand side is itself an integer. Thus, SJ 1 C V fc n and: 



p(U^r) < p(U ^ n ) • ( 5 - 3 ) 

We define also the sets: 

Let u e {V? \ WJt). Then by h^^,^ < h^ )1/d]+ we get: 

i : ^ [((fe+1)can7)V ^ < M> ttli : ^' [(feCO " 7)1/d]+ < - > > k + 1 . 



n ■> i- n 
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Hence, (V fc n \ W£) C V fc n | _ 1 , and therefore we have for any fixed n and k: 

V k n C W£ U V k n +1 . (5.5) 
Applying this inclusion M times, for k — 1, . . . , M , we obtain: 



M M M / M 



U^Tc (Vf u|Jn n ) c (H/fu^u|J^ n ) c-.-c |J^ n uK m+i • ( 5 - 6 ) 

k=l k=2 k=2 \k=l J 

Then we take the limit M — > oo to recover the infinite union that one needs in (15.31) and we 
use the inclusion ( 15.61) to find the inequality: 

M 



fe>l fc=l 

M oo 

< lim (^P(^) + P(^ +1) )) = J2H W k) + I/™ P(VS). 



M^oo v 1 T ' / v 7 M^oo 

fc=l fc=l 



The limit in the last term can be calculated directly : 



J+ <- \>M) (5. 



lim P(V£) = lim pJ w : ti/i : # 
= lim P< u : Qn7u/dl (1/n) ^ — T \ = F{ u : f(l/n) ^ AV 

A/^oo 1 [(A/e-" 7 )i/d] + V / y ^ [(Me Qn7 ) 1 /^] rf 

for some constant K. In the last step we used dominated convergence theorem. 

Now we can use the Lifshitz tails representation for the asymptotics of the a.s. non-random 
limiting IDS, v(E), see ( 13.111) . which implies: 

lim sup e and/2 v{l/n) < 1 , (5.9) 

n— >oo 

for a > 0. Since we assumed that < 7 < d/2, there exists n < 00 such that by (I5.8P and 
(15.91) for all n > no we get: 

lim ¥{Vm) = 0. 

M— >oo 

This last result, along with ( 15. 3ft and ( 15.71) . implies that: 

00 

E P (U 5 < E E p (^ n )- ( 5 - 10 ) 

n>no Z>io n>no fc=l 

Now, we show that the upper bound in (15.101) is finite. First we split the box A^ feea „ 7 )i/dj 
into m(k,n) disjoints sub-cubes of the side l(k,n), with the following choice of parameters: 

m(k,n) := [kM n } + , M n := B~ d/2 e anl n~ d/2 
[{ke an ~') l ' d ] 



l(k, n) 



(m(k,n)y/ d 
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Here B is the constant that comes from Lemma 15.11 Now by the Dirichlet-Neumann in- 
equality, see e.g. [T3], we get: 

m(k,n) 
+ 3=1 

where hf^ n ^ denotes the Schrddinger operator defined in the j-ih sub-cube of the side l(k, n), 
with Neumann boundary conditions. Note that, by the positivity of the random potential, 
we obtain: 

E if > 4 > 777^ > ~ ■ (5-12) 



Z(fc, n) 2 n 

l{k,n)i 



Here E^' 2 denotes the second eigenvalue of the operator h^, kn y and e^ 2 the second eigenvalue 

of — Aj/^x , i.e. the kinetic-energy operator defined in the j-th sub-cube of the side l(k,n) 
with the Neumann boundary conditions. 

By equation (I5.12p . we know that to estimate the probability of the set (j5.4p by using 
the Dirichlet-Neumann inequality (15.111) . only the ground state of each operator hf^ n ^ is 
relevant. Since the sub-cubes are stochastically independent, we have: 

P(W£) < pjo; : ${j : E%f < 1/n} = A; J < ^^/(l-g)*'^ < m ^C k q k 

with g being the probability P{c<j : < 1/n}. The latter can be estimated by Lemma 

15.11 So, finally we obtain the upper bound: 

F(W£) < exp{-/L4(Z(fc,n)) d } . (5.13) 

Using Stirling's inequalities, see [T6] : 

(27r) 1 / 2 n n+1 / 2 e- n < n! < 2(27r) 1 / 2 n n+1 / 2 e - n . 

we can give an upper bound for the binomial coefficients m ^ k ' n ^Ck in the form: 

2{2n)^(kM n + 5)(kM n +6+i/2) exp (_fc Mn + 5) 
(27r)fc fc+ 5 exp(-fc) • (kM n + 5- fc)(fcAf„+<s-fc+i/2) eX p(-A;M n + 5 - Jfe) ' 

where 5 ^ is defined by: 

m(A;,n) = [A;M n ] + = kM n + 5 . 

Then (15.141) implies the estimate: 

„«,„ )Ct s (M<, + f)^V» , (1 + ,)^, 

k k+ ^(kM n - k) kM ^ & - k+1 / 2 V(l - (72 )( fcM «+ l5 +3- i: ) 

for some i^i > and 

^ := <5(A;M n ) _1 , a 2 := M" 1 . 



(5.14) 
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Since 8/k < 1 and cr^ — > as n — > oo, and also using the fact that xln(l + — ► 1 as 
x — > oo, we can find a constant c > such that, for n large enough one gets the estimate: 

m(M) ^ < ^( M ") fc ( (l- M-")(^-^) ) ^ A ^M„)V fc . (5.15) 
The side l(k,n) of sub-cubes has a lower bound : 

l(k n) ~ [(fceQ " 7)1/d] + > ^ 1/d > (B«* -±-\ 1M (5 16) 
Combining (I5.15p . (15.161) and (15.131) we obtain a sufficient upper bound: 

< K 2 J^exp jfc(cm 7 - (tZ/2) \n(nB) + c - AB d/2 n i,2 \ } 



-fcAZ d (fc,n) 
\ fc a ck a -kAB d / 2 n d / 2 /(l+o-i) 

k>l 



k>l 



^ K 3 J2 ex P *: («™ 7 - ^5 d/2 n d/2 + if 4 ) ^ K 5 exp{-K 6 n d/2 ) . 



Here Ki are some finite, positive constants independent of k, n, I, for any n large enough. 
Now the lemma immediately follows from (15.101) . □ 

Proof of Theorem 15.11 Let A n to be the event: 

A n := {u : if (1/n) > e~ an ~' for some /} . (5.17) 

By Lemma [5.21 we have: 

£>(4») < oo, 

and therefore, by the Borel-Cantelli lemma one gets that with probability one, only a finite 
number of events A n occur. In other words, there is a subset Q C Q of full measure, P(fi) = 1, 
such that for any uj G Q one can find a finite and independent on Z number no(c^) < oo for 
which, in contrast to (15.171) . we have: 

if (l/ra) e" Qn7 , for all n > n (uo) and for all Z ^ 1. 

Define := l/no(uj). For any ^ S(u), we can find n ^ no(o;) such that: 

— ^ E ^ -, 
2n n 

and the theorem follows with the constant a modified by a factor 2~ 7 . □ 
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6 On the nature of the generalized condensates in the 
Luttinger-Sy model 



In this section, we study the van den Berg-Lewis-Pule classification of generalized BE con- 
densation (see discussion in Section |3} in a particular case of the so-called Luttinger-Sy 
model with point impurities [3J. 

Let u(x) > 0, x G R, be a continuous function with a compact support called a (repulsive) 
single-impurity potential. Let {f^x}^^ be the random Poisson measure on R with intensity 
A > : 

F({w6fl: ^(A) = n}) = e" A|A| , n e N = N U {0} , (6.1) 

n\ 

for any bounded Borel set Acl. Then the non-negative random potential generated by 
the Poisson distributed local impurities has realizations 



v (x\ 



j ^(dy)u(x-y) = ( 6 - 2 ) 



Here the random set X w corresponds to impurity positions X w = {^"} C R, which are 
the atoms of the random point Poisson measure, i.e., Jj {X u \ A} = ^(A) is the number of 
impurities in the set A. Since the expectation E(z/^(A)) = A |A|, the parameter A coincides 
with the density of impurities on R. 

Luttinger and Sy defined their model by restriction of the single-impurity potential to 
the case of point 5-potential with amplitude a — > +oo. Then the corresponding random 
potential (16.21) takes the form: 

I v»(dy)a8(x-y) = a ^ 5(x - a%) . (6.3) 

Jr ijgx" 

Now the self-adjoint one-particle random Schrodinger operator h% := h° + v^ is defined in the 
sense of the sum of quadratic forms (12. 2\i . The strong resolvent limit h°[ s := s.r. lim a ^ +00 
is the Luttinger-Sy model. 

Since X^ generates a set of intervals {/" := (x^_ x , x^ )}. of lengths {L^ := Xj — sc^}., 
one gets decompositions of the one-particle Luttinger-Sy Hamiltonian: 

^L5 = 0M^). dom(^)C©I 2 (/J) , WGO, (6.4) 



into random disjoint free Schrodinger operators {^dC^/)} with Dirichlet boundary con- 
ditions at the end-points of intervals {-^} • • Then the Dirichlet restriction hf D of the 
Hamiltonian h^ s to a fixed interval A; = (—1/2,1/2) and the corresponding change of nota- 
tions are evident: e.g., {/?}. ^ {lf} M = ^\ where M' (a>) is total number of subintervals in 
A; corresponding to the set X^ . For rigorous definitions and some results concerning this 
model we refer the reader to [S]. 

Since this particular choice of random potential is able to produce Lifhsitz tails in the 
sense of (13. lip , see Proposition 3.2 in [5], it follows that such a model exhibits a generalized 
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BEC in random eigenstates, see (13.41) . In fact, it was shown in [5] that only the random 
ground state <^' z of hf D is macro scopically occupied. In our notations this means that 

lim = { - [t f lt Pc (6-5) 

lim - = 0, for alH > 1 . 

I— >oo I 

According to the van den Berg-Lewis-Pule classification this corresponds to the type I Bose- 
condensation in the random eigenstates {^ji^i. 

Following the line of reasoning of Section HI we now consider the corresponding BEC in the 
kinetic-energy eigenstates. We retain the notation used in that section and explain briefly 
the minor changes required in the application of our method to the Luttinger-Sy model. 

We first state the equivalent of Theorem 14.11 for this particular model. 



Theorem 6.1 Theorem\4- 1\ holds with the function g defined as follows 



e -||x|| 2 (l/2n/3) 



(2nn(3) d / 2 

/ w n *{dt) exp ( - A(sup£(s) - inf £(*))). 



x 

' So,*) 



The scheme of the proof is the same as above, cf. Sections H] and El First, we note that 
Lemmas 14.11 and 14.21 apply immediately. The positivity of the random potential has to be 
understood in terms of quadratic forms, see (12. 4p . 

Before continuing, we need to highlight a minor change concerning the finite-volume Lif- 
shitz tails arguments. Although the Theorem 15.11 is valid for the Luttinger-Sy model, its 
proof (see Section [3]) requires a minor modification, as the assumption of Lemma 15.11 is 
clearly not satisfied for the case of singular potentials. However, by direct calculation we can 
obtain the same estimate with the constant B = 7r 2 /4 in (15.11) . First, suppose that there is 
at least one impurity in the box, then the eigenvalues will be of the form (for some j) 

(nV)/(LJ) 2 , 71=1,2,... 

if Vj is an inner interval (that is, its two endpoints correspond to impurities), and 

((n + l/2) 2 7r 2 )/(LJ) 2 , n = 0,1,2,.... 

if If is an outer interval (that is, one endpoint corresponds to an impurity, and the other 
one to the boundary of A;). Therefore, E"' l,N ^ B/l 2 since obviously < /. Now, if there 
is no impurity in the box A/, then Ef' ' N = < B/l 2 . But due to the Poisson distribution 
(16. ip this happens with probability e~ xl , proving the same estimate as in Lemma [5.11 
With this last observation, the proof of the Theorem 15.11 in Section [5] can be carried out 
verbatim, without any further changes. 

Our next step is to split the measure rhi into two, and m[ 2 \ see (I4.9p . and prove the 
statement equivalent to the Theorem 14.21 
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Theorem 6.2 For any d ^ 1, the sequence of Laplace transforms of the measures 



/,(*;/?, := / mf^de) e" te 
converges for any t > to a (non-random) limit f(t; fiiHoo) > which is given by: 

__ /• -\\x\\ 2 (l/2n/3+l/2t) 

(df) 8xp(-A(sup{(»)-inf {(»))). 



a/ 

(O,a0 



X / W 



Proof: We follow the proof of Theorem 14.21 using the same notation. The uniform conver- 
gence is obtained the same way, since the bounds (j4.14j) . (14.151) . and (14.161) are also valid in 
this case. As in (14.201) . we can use the ergodic theorem to obtain: 



r -\\x\\ 2 (l/2np+l/2t) r 

lima,(n) =e*»E„ 6* /2 g «^(d£) WO- («*) 

j (0,a:) 

We have used the fact that the Dirichlet boundary conditions at the impurities split up the 
space Tii into a direct sum of Hilbert spaces (see (16.41) ). This can be seen from the expression 

.||" (l/2„/9 + l/2t) , -rd S aE,. e x^«( S )-^) 



/" p -||x|r(l/2n/^+l/2t) /• 

lim a,(n) = e""- / dx E w / (d£) e 



by formally putting the amplitude, a, of the point impurities (I6.3P equal to +00. Because of 
the characteristic functions Xi^,n/3i which constrain the paths £ to remain in the interval 
in time n/3, the sum in (16. 6p reduces to only one term: 

p — ||x|| 2 (l/2n/3+l/2i) r 

fan = tf*- / dx 6 J/2 IS. y x (0 „, M ,„Mf), (6.7) 

(0,1) 

where (a^, is the interval among the J"'s which contains 0. 

The expression in (16.71) can be simplified further by computing the expectation explic- 
itly. 

First, note that the Poisson impurity positions: a w , b u are independent random variables and 
by definition, a u is negative while b^ is positive. For the random variable b^ the distribution 
function is: 

IP {b u < b) := P{(0, b) contains at least one impurity} = 1 — e~ Xb , 
and therefore its probability density is Xe~ Xb on (0, 00). Similarly for a w one gets: 
P (a w < a) := P{(a, 0) contains no impurities} = e" A|a| = e Aa , 
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and thus its density is \e Xa on (— oo,0). Using these distributions in (16.71) we obtain: 

/0 poo p -||x|| 2 (l/2n/3+l/2t) 

dae Aa / dbe~ Xb / dx — x 

-oo J0 JR 



(A7rHnf3) d / 2 



fO roo /• -\\x\\ 2 (l/2nf3+l/2t) 

; n ^°°A 2 / dae Aa / dbe~ Xb / dx - - x 



x / w n \di) l(sup(£(s)) ^ b) l(inf(£( S )) ^ a) 



-\\x\\ 2 (l/2nfl+l/2t) 

3 n ^°°A 2 / dx———^—x 
(AiT 2 tn(3) d l 2 



X 



dae Aa / d6e- Ab 



and the Theorem 16.21 follows by explicit computation of the last two integrals. □ 

Proof of Theorem I6.lt Having proved Theorem 16.21 it is now straightforward to derive 
the analogue of Corollary 14.11 for the Luttinger-Sy model. Note also that the Corollary 14.21 
remains unchanged, since only the uniform convergence was used. With these results, the 
proof of Theorem 16.11 follows in the same way as for Theorem 14.11 □ 



We have proved, in Theorem l6.ll that the Luttinger-Sy model exhibits g-BEC in the kinetic 
energy states. But, in this particular case, we can go further and determine the particular 
type of g-BEC in the kinetic energy states. Recall that the g-BEC in the random eigenstates 
is only in the ground state, that is, of the type I, see (16.51) and [5] for a comprehensive review. 
Here we shall show that the g-BEC in the kinetic-energy eigenstates is in fact of the type 
III, namely: 



Theorem 6.3 In the Luttinger-Sy model none of the kinetic- energy eigenstates is macro- 
scopically occupied: 

lim \ (Ni(il) k ))i = for all k E A?, 

I— >oo I 

even though for ~p > p c there is a generalized BEC. 



To prove this theorem we shall exploit the finite- volume localization properties of the random 
eigenfunctions of the Hamiltonian hf D . Since the impurities split up the box A; into a 

finite number M l (cu) of sub-intervals {if}*^^, by virtue of the corresponding orthogonal 

decomposition of hf D , cf (16.41) . the normalized random eigenfunctions are in fact sine- 
waves with supports in each of these sub-intervals and thus satisfy: 

W\x)\< , l<j s <M\u). (6.8) 

We require an estimate of the size Lf of these random sub-intervals, which we obtain in the 
following lemma. 
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Lemma 6.1 Let A > be a mean concentration of the point Poisson impurities on R. Then 
eigenf unctions (p'j are localized in sub-intervals of logarithmic size, in the sense that for any 
k > 4, one has a.s. the estimate: 

hmsup ^-p— --— ■ — - ^ -. 

i^oo In I A 



Proof: Define the set 



Si := {uj : max > — In/). 



Let n := [2A//(/«ln/)] , and define a new box: 

Split this bigger box into n identical disjoints intervals {I l m }m=i °f s i ze ^(2A) _1 ln/. If 
uj G Si, then there exists at least one empty interval I l m (interval without any impurities), 
and therefore the set 

Si C [J {uj : I l m is empty} . 

By the Poisson distribution ( 16. II) . the probability for the interval to be empty depends 
only on its size, and thus 



) ^ n exp(-A— In/) < 



2A/ 
Kin/ 



Since we choose k > 4, it follows that 



i) < oo. 



Therefore, by the Borel-Cantelli lemma, there exists a subset Q C f2 of full measure, P(fi) 
1, such that for each a; G f2 one can find Iq(uj) < oo with 



P {uj : max Z# < - In/} = 1 



A 



for all / ^ /o(<^) 



□ 



Now we can prove the main statement of this section. 
Proof of Theorem 16. 3t The atom of the measure fh has already been established in 
Theorem 16.11 Concerning the macroscopic occupation of a single state, we have 



dx ip k {x)$' l (x) 



7 £<W')>i 



x)| 
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where in the last step we have used the bound ^ \j\fl . Therefore, by (16. 8ft and Lemma 
16.1^ we obtain a.s. the following estimate: 

i 

which is valid for for large enough I and for any k > 4. The theorem then follows by taking 
the thermodynamic limit. □ 



7 Application to weak (scaled) non-random potentials 

It is known for a long time, see e.g. p2], [IE], that BEC can be enhanced in low-dimensional 
systems by imposing a weak (scaled) external potential. Recently this was a subject of a 
new approach based on the Random Boson Point Field method [19]. In this section, we show 
that, with some minor modifications our method can be extended to cover also the case of 
these scaled non-random potentials. 

Let v be a non-negative, continuous real-valued function defined on the closed unit cube 
Ai C M. d . The one-particle Schrodinger operator with a weak (scaled) external potential in 
a box A/ is define by: 



2 



A D + v(xt/l, . . . , x d /l) . {1.1) 



Let {(p\, E\}i^i be the set of orthonormal eigenvectors and corresponding eigenvalues of 
the operator (17.11) . As usually we put E\ ^ E2 ^ . . . by convention. The many-body 
Hamiltonian for the perfect Bose gas is defined in the same way as in Section [21 We keep 
the notations m and m for the occupation measures of the eigenstates and of the 

kinetic-energy states respectively. We denote the integrated density of states (IDS) of the 
Schrodinger operator (17.11) by ui, and by v = lim^oo v\ its weak limit. We assume that the 
first eigenvalue E[ — > as I — > 00, which is the case, when e.g. v(0) = 0. This assumption 
is equivalent to condition (iii), Section [2j It ensures that for a given mean particle density p 
the chemical potential ii OQ {(5 ) 'p) satisfies the relation (13.61) . where Ji := 7i(/9,p) is a (unique) 
solution of the equation [T7] : 

for p < p c , where the boson critical density is given by: 

Here u° is the IDS (13.91) of the kinetic-energy operator (12.11) . In particular the value p c = 00 
is allowed in (17.3j) . If p c < 00, the existence of a generalized BEC in the states {(p l j}i>i 
follows by the same arguments as in Section 3. For example, the choice: v(x) = \x\, makes 
the critical density finite even in dimension one, see e.g. [T7] . 

Now, we prove the statements equivalent to the Theorem 14.11 
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Theorem 7.1 The sequence {rhi}i^i of the one-particle kinetic states occupation measures 
has a weak limit m given by: 



C ) 



m[Ci£) ~ \ F(ey(de) , ifp<p c , 

where the density F(e) is defined by: 

F(e) = [ dx (e^+<*)-»°°) - i)- 1 ; 
and /Zoo := Poo(fl,p) satisfies the relation ( Iff. 6]) . 

We note the similarity of this result with the free Bose gas. Indeed, the kinetic-energy states 
occupation measure density is reduced to the free gas one, with the energy shifted by the 
external potential v and then averaged over the unit cube. 

The proof requires the same tools as in the random case. As before, we split the occupation 
measure into two parts: 

m = fh^+mP with 
mf\A) := VihP^e-^jl^Kl/n), 

rhf\A) := Vi T rP A ( e -*^)l( tt >l/n), 
and we prove the following statement: 

Theorem 7.2 The sequence of measures fn^ 1 converges weakly to a measure m^ l \ which is 
absolutely continuous with respect to z/° with density F(e) given by: 



F(e) = [ dx ( e /^+^)-M~) _ i) 1 



Proof: We follow the line of reasoning of the proof of Theorem 14. 21 Let g\(t\ /3,pi) be the 
Laplace transform of the measure mf 1 -. 



gi{t;P,m) = / mf\dE) e" te (7.4) 
Jr 

= ^ Tr e"^ ? { e - nP{hl ~^) l(/x, ^ l/n) 



Again, our aim is to show the uniform convergence of the sum over n with respect to I. Let 

^-Tre~ th ° e- n/3(/l '- w) l(/z z ^ l/n) (7.5) 



ai[n) :-- 

Vi 



^ ^Tre- n/3{hi - w) l(/i Z ^ l/n) 
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Then for p < p c we can apply a similar argument as for the random case, since the estimate 
\i>i < /Woo/2 < still holds, to obtain: 



ai(n) ^ e n/W2 I e^v^de) ^ 1^ e «/W2 

J[0,oo) 



[0,oo) 

If p ^ p c , then pLi ^ 1/n in (17.51) implies that: 

a*(n) ^ e^V e" n ^ ^ - — ^— — / dxe~ n ^ x \ 
Y (2vrn/3) d / 2 i Al 

where the last estimate can be found in [T7] or [IB]. Now the uniform convergence for the 
sequence ai{n) follows from (17.31) . since we assumed that p c < 00. The latter implies also 
that for p > p c , /ioo(/5, p) = 0. Thus, we can take the limit of the Laplace transform (17.41) 
term by term, that is: 

lim a z (n) = lim — Tr e~ th ° e- nl3{h '-^ l) l{pi ^ 1/n) (7.6) 

l— >oo I— >oo Vi 

= lim 1 / / dxdx'e^^^xOe-^-w^x',^) 

Iff e -\\x-x'\\ 2 (l/2n(3+l/2t) 

e»*»lim^ / / dxdx' x 



^ A, A, (4vr 2 tn/^/ 2 



nil 



(ic,:r' J (x' ,x) 

th?i 



Here we have used the Feynman-Kac representation for free e - * i(x, y) and for non-free 
e -f3h t ^ Qj^^g semi-group kernels, where w T stands for the normalized Wiener measure 
on the path-space Qj x y s , see Section 14.11 

Note that by Lemma IA.21 which demands only the non-negativity of the potential v, we 
obtain for (17. 6p the representation: 

lim —Tre- th ° e ~ n ^-^) (7.7) 



1^00 Vi 



r 1 P „-\\x-x'\\ 2 (l/2n/3+l/2t) [■ 

V ^ 7a, (4vr 2 tn/3) d / 2 V, _ V W 



(4vr 2 tn,5) d / 2 

(x f ,x) 



Now we express the trajectories £ in terms of Brownian bridges «(t) e f2, ^ r ^ 1, we 
denote the corresponding measure by D. Letting x = x'/l, we obtain: 

lim —T Te - th i e- n(3{h '~^ 



V, 



e -\\x-lx\\ 2 (l/2nl3+l/2t) 
D nt3tMca lim / dx / dx t— — - 777: X 



V Ax (4vr 2 tn/3)<V 2 

x [ D(da)ex V (- T dsv[(l - ^-)x + -^(x/l) + ^-a(s/n(3)}) . 
in ^ Jo n P n P ' ' 
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Since the integration with respect to x is now over the whole space, we let y = x — Ix to get 
lim — Tr e~ th ° e~ n ^ hl ~^ 

l— >oo vi 



JRd J Al 



e -||y|| 2 (l/2n/3+l/2t) 



(A7T 2 tn/3) d / 2 

n(3 



X 



x / D(da) exp ( — / dsv[x-\ ^(y/l)-\ -, — a(s/n(3) 

Jn ^ Jo n P ' 

f -\\y\\ 2 (l/2n(3+l/2t) f 

J Rd (Stripy/ 2 J Al 

where the last step follows from dominated convergence. Therefore, we obtain by (17.41) the 
following expression for the limiting Laplace transform: 

inn 9l \t, p, m) 2^ e (2vr ( n/9 + t ))d/2 J axe 



It is now straightforward to invert this Laplace transform (for each term of the sum), to find 
that: 

F(E) u\dE) = lim rh\{dE) = ^V"^ - " 00 * ( f dxe~ npv{x A u°(dE) . 

The Theorem then follows by Fubini's theorem. □ 

Proof of Theorem I7.lt The proof of Theorem 14.11 can be applied directly. Note that 
Lemmas 14.14 14.21 are still valid, since (as we emphasized in Remark 14. lj) their proofs require 
only the non-negativity of the external potential. Similarly, Corollary 14.21 now can be used 
directly, since we have proved Theorem 17.21 □ 



Appendices 

A Brownian paths 

In this section, we first give an upper estimate of the probability of a Brownian path to leave 
some spatial domain, cf. e.g. [20] and the references quoted therein. 



Lemma A.l Let the set 

nf^o : = = £(0) = x, £(T) = x'} 

be continuous trajectories from x to x' with the proper time < r < T , and with the 
normalized Wiener measure w T on it. Let x,x' be in A/, and Xa,,t(£) the characteristic 
function over OT,^ of trajectories £ staying in A/ for all ^ r ^ T. Then one gets the 
estimate: 

I ^(d£)(l- XAi , T (£)) ^ e -c ( T)( m inM^),^'^)}) 2 . (A1) 

(x,x') 
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Proof: Define a Brownian bridge a(s),0 ^ s ^ 1 by: 

f (t) = (1 - r/T) x + t/T x' + VTa(r/T). 

Let us consider first the one dimensional case, i.e. A; = [—1/2, 1/2]. Without loss of generality, 
we can assume that: 

d(x,dAi) «C d(x',dAi). 

Suppose that x > 0, then we have: 

—x ^ x' ^ x and <9A;) = 1/2 — x 

Assume that the path £ leaves the box on the right-hand side. Then, for some t, we have: 

m > l - 

<*(t/T) > ^(^ + (t/T-l)x-^ 

ait/T) > l(L + ( t /T-l)z-±z)=±d(x,dA l ) (A.2) 



/TV2 w ' T ) y/f 

The case, when £ leaves the box on the left-hand side can be treated similarly. 

Let x < 0, then we have: 

x ^ x' ^ —x and <9A Z ) — 1/2 + x 

Again, assume that the path leaves the box on the right hand-side. Then, for some t, we 
have: 

m > \ 

<*(t/T) > -j=(t + ( t /T-l)x-±;J 

a(t/T) > -L(L-(t/T-l)x'-j;x') > -j=d{x,dh) (A.3) 

The case, when £ leaves the box on the left hand-side can be considered similarly. The 
relations flA.2j) . (IA.3[) imply that if £ leaves the box A; in one dimension, then the Brownian 



bridge a must satisfy the inequality: 

sup | a (t/T) | > C(T)mm{d(x,dAi),d(x',dAi)}, (A.4) 
t 

for some constant C(T). 

This observation can easily be extended to higher dimensions, when x := (xi, . . . , xj) and 
a(s) := (cti(s), . . . , arf(s)). Now, if £ leaves the (ti-dimensional) box A;, there exists at least 
one i such that similar to (IA.4I) : 

sup | oa (t/T) | > C(T)mm{d(x l ,d i A l ),d(x' i ,d l A l )}, 
t 
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where we denote d(xi,diA{) := min{//2 — Xi,l/2 + Xi}. Now, since A/ are cubes, we get 
d(xi, diA[ ^ d(x, dAi for any x G A;. Then we obtain: 

\\a(t/T)\\ > \ai(t/T)\, i = l,...,d, 

sup ||a(t/T)|| > maxsup \cti(t/T)\ , 
t » t 

sup||«(t/T)|| > C^min^,^),^,^)} 
t 

> C(T) min{d(x, 9A,), d(x', OAfi} . (A.5) 

Therefore, the probability for the path £ to leave the box is dominated by the probability for 
the one-dimensional Brownian bridge a to satisfy (1A.5D . The latter we can estimate using 
the following result from [20] : 

p(supa(s) > x) ^ Ae~ Cx2 
valid for some positive constants A, C, which implies the bound ( lA.ip . □ 



Now we establish a result, that we use in the proof of Theorem 14.21 

Lemma A. 2 Let K^^x^x'), Kqi(x,x') } Kq(x,x') be the kernels of operators exp(-thf), 
exp(-th^), and exp(— t A/2) respectively. Then 

lim \- I [ dxdx , K% l {x,xf)Kf l {J,x) (A.6) 

l^oo Vi J Al J A[ 

= hm [ dx±- [ dx'K^&x') [ w n ^dOe~^^ {i (s))_ 

l->oo J Rd Vl J A J n f 

(x f ,x) 



Proof: By the Feynman-Kac representation, we obtain: 



lim — 

l— too Vl 



l - I I dxdx'Kl^x')K n ^{x\x) 

V l J A, J A, 
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1 r r „-\\x-x'\\ 2 (l/2n(3+l/2t) 



Vi J a, J ^ (4nHnpy/i J Q ne x) 

in* ,. 

To eliminate the characteristic functions restricting the paths in the last integral, we 
shall use Lemma [A. II First, we estimate the error 7(d) when we remove the restriction on 
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the path £: 



If f e -\\x-x'\\ 2 (l/2n/3+l/2t) 

7(d) := lim — dx dx' t— — — To x (A. 7) 



/ 

n0 



w^(dO e-f° ds ^» (1 - XmMO) / w*(dO WO 

1 /" r -||x-x'|| 2 (l/2n/3+l/2t) /• 

1 /■ /■ e -||x-x'|| 2 (l/2n/3+l/2t) 

« , h i h A d * A ^w^ A <} > *'. SA ')> (4^)^ x 



ti^dfl (1 - xWO) 

I /■ /■ -\\x-x'\\ 2 (l/2nf3+l/2t) 



1 /" /" e -||x-x'|r(l/2n^+l/2tj 

m »"(df) (1 - X A„nMe) 



n(3 
(3:', 3:) 



^ lim 1 / dx / dx / ^(x,x , )/C /3 (^^)e" C(n/3)(d{x '' 9Ai)2 

/ — s-OO V/ 



A, «/A 



+ lim ^- f dx [ dx'Kl{x,x')Kf{x',x)e 
1^00 Vi J A J A 



■C(n/3)(d(x,dAi) 2 



where the last step is due to Lemma [A. 11 Since all integrands are positive, we can extend 
one of the spatial integrations to the whole space, and hence we get: 



7(d) ^ lim \- [ dx [ dx' Kl{x,x')K { 

l-*oo Vi J R d J A 



+ lim 1 / dx / dx'K t (x,x')K^(x',x)e- c ^ )Wx ' dAl)2 
1-+°° Vi J Al J R d 

^ lim -Kl +nP [ dx > e -c(nm(x'Mi) 2 
Vi J Al 

+ lim l-Kl +nP [ dxe-^W^', 

l^oo V t J Al 

where we have used the notation i\~Q +n/3 := Kq +uI3 (x,x) since these are independent of x. 
Finally, using the fact that the boxes A; are cubes of side I, we obtain: 



T (d) <C lim ^!L_ / dx ' e -C(nm/2-x'f + Um f^0_ / dxe -C(n/3)m-xf = Q 

I-+00 I J-1/2 I J-l/2 
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We can estimate the error estimate due to the removal of the characteristic function for £' 
in (I4.18P in the same way. Therefore, we get: 

Iff e -\\x-x'\\ 2 {l/2nP+l/2t) 

lim — / / dxdx' ——z — —jjz x (A. 8) 

f-n/3 



Jnf, . Ju* 

VX ,X) \x,x' J 

1 r r p -\\x-x'\\ 2 {l/2nf3+l/2t) r p 

lim- / / dxdx' , „ — m / w ^(dOe-^ ds ^^) uj^fdf) 



Now we show that one can replace the first integration over the box A; by one over the 
whole space. Let 7(d) be the error caused by this substitution. Then by the positivity of 
the random potential we get the estimate: 



If f e -\\x-x'\\ 2 (l/2n(3+l/2t) 



X 



cif, , J si* 

[x 1 ,x) {x,x' ) 



If f e -\\x-x'\\ 2 {l/2nl3+l/2t) 

< lim — / dx / dx' 



1^00 V, J Rd ^ Ai J Ai (AnHnpyi* 
In the one-dimensional case the estimate of the error term (1A.9j) takes the form: 

1 /— «/2 fl/2-x -y 2 (l/2n/3+l/2t) 

,(1) * £l n 7 ^ ds/^d, (4 ^ )1/2 (A,0, 

/-oo /.J/2-X e -j/ 2 (l/2n/3+l/2t) 

+ lim — / dx dy 



* A 2 V-i/2-x (4vr 2 tn/3) 1 /2 
For the first term one gets: 

I f-l/2 fl/2-x e -y 2 {l/2nf3+l/2t) 

lim — / dx / dy 



I Aoo y-i/2-x (47r 2 tn/5)V2 

1 f l e -y 2 (l/2nP+l/2t) rl/2 j /-oo e -j/ 2 (l/2n/3+l/2i) W/2-J/ 

lim - / dy — 7— - — -n / dx + lim - / dy — 7— - — ^ , n / dx = . 

One obtains a similar identity for the second-term in ( 1A.10I) . Direct calculation shows that, 
the error term for higher dimensions ( 1A.9I) reduces to a product of one- dimensional terms 
(fATTOl) . Then ([Oil gives: 

Iff e -||x-x'|| 2 (l/2n/3+l/2t) 

lim — / / dxdx' -— — -—777; x (A.ll) 

Vi J Al J Al (Stripy/? v ; 



x / w^(dO *•«"<««) ^ (£) / ^(d^O XAj , t (0 

(ar ,ac; {x,x' ) 

r 1 r „-\\x-x'\\ 2 (l/2nf3+l/2t) r . 

lim / dx-W dx' 6 j— / ^(dOe-^-^)) 



which finishes the proof of (IA.60 . □ 
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B Some probabilistic estimates 



First we recall the assumptions on the random potential v w used in [9], and which we also 
adopt in this paper: 

1. (a) On the probability space (Q, T , P) there exist a group of measure-preserving met- 

rically transitive transformations {T p } peR d of Q, such that v u (x+p) = v TpUJ (x) for 
all x,p E R d ; 

(b) Ec;{/ Al dx \v"(x)\ K } < oo, where k > max(2, d/2). 

2. For any A C M d , let Sa be the cr-algebra generated by the random field v UJ (x),x G A. 
For any two arbitrary random variables on Q, f,g satisfying (i) \g\oc < oo, E w {|/|} < oo 
and (ii) the function g is E^i -measurable, the function / is EA 2 -measurable, where 
Ai, A2 are disjoint bounded subsets of M d , the following holds 

\E{\f.g\} - E{|/|}E{|<7|}| < |^| oo E{|/|}0(d(A 1) A 2 )) 

with (j){x) — > as x — > 00, and rf(A 1; A 2 ) the Euclidean distance between Aj and A 2 . 

After recalling these conditions, we can give a sketch of the proof of Lemma 15.11 

Let h^' N to be the Schrodinger operator ( 12. 3ft . with Neumann boundary conditions instead 



of Dirichlet, and denote by {E"' ' its ordered eigenvalues (including degeneracy) 

and the corresponding eigenvectors. Similarly we define the kinetic energy operator fi}' N with 
the same boundary condition, and denote by {e l £ , ip l ^ N }k^i its ordered eigenvalues (including 
degeneracy) and corresponding eigenvectors. The following result is due to Thirring, see [21] : 



Lemma B.l Let v^ a := v u + Xa, for X,a > 0. Then, 

E»M > -\ a + minj^, 1 jf dx KJx))' 1 } . 

Proof: Let P to be an orthogonal projection in J^. Then for any vector from the 
intersection Q(v^ a ) f]Q((v^ a ) 1/2 P(v^ a ) 1/2 ), we have: 

(0X,» = («,J 1/2 0,«J 1/2 0) 

= («J V V, ^KJ 1/2 0) + («J 1/2 0, (i - ^)KJ 1/2 0) 
> (KJ 1/2 ^KJ 1/2 0), 

and therefore, 

- ^ + < a £ -£A* + KJ 1/2 P«J 1/2 , (B.l) 
in the quadratic-form sense. Let us choose: 

p ■= (^-^((^.KJ-vr))" 1 ^-)- 172 
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where P is the orthogonal projection onto the subspace spanned by the vector . It can 
be easily checked that P is an orthogonal projection. Applying fIB.ip to the function <j)^ l,N 
one gets: 



Ef' N + Xa > 



l lo,1,N ( i A \iuj,1,N\ . \/,lu,1,N i l : N \ 1 2 ( il,N / w \-l /iJV\— 1 



'1 ) 



s \/iw,l,N il,N\\2 l,N . \/iw,l,N 1 1, 

> 2^ \m )l £ k + A 



JV\i2 



fe>l 



1 



dsKJ*))" 1 



-i 



But since the Neumann boundary conditions imply that e{ N = 0, we obtain 



E x " + Xa ^ (1 - K^" ,ip{ )| )e 2 ' +1(0!" ,t/j 1 



u>,l,N I l,N\ 1 2\ Z,AT 



1 



dxK-Cx))- 1 



To finish the proof, we have to study separately the two cases, namely, e 2 ' l ess than and 

-l 

□ 



greater than 



Proof of Lemma 15.11 By Lemma [B. ll with A = B/Z 2 and a as defined in assumptions, i.e. 
for B — 7r/(l + a), a > p/ (1 — p), we have: 



Ef' N > 



aB 
l 2 



+ min(7r// 2 , 



where X 1 ! 



dx 



1 



Vi J Al + 5a// 2 



Therefore, 



£ -^ + min(vr// 2 ,l/Xr). 



5 



7T 



Hence, the inequality E%' 1 ' < B/l 2 implies that Xf > 1 2 /tt and consequently: 



F(E^ l ' N < B/l' 2 ) < F(X? > 1 2 /tt) 



(B.2) 



Define a random variable ^(5) := Vi 1 J. dx 6/(v w (x) + 5), which is an increasing function 
of 5. Then for the left-hand side of ( IB. 21) one gets the estimate: 



¥(Ef' N < B/l 2 ) < P(V(Ea/Z 2 ) > - 



a 



By Lemma 2 in [9], we know that for any positive S the random variables {Y^(5)}z, converges 
geometrically to a limit 1^(5) as / — > oo, that is, for any e > 0, there exists a constant M(S, e) 
such that 



\Yf(6)-Y 00 (6)\>e/2) < e 



-M(<5,e) V5 



for / sufficiently large. By the ergodic theorem Y^S) is non-random and can be expressed 
as: 



Y^d) = E u 



5 

v"(0) + 5 
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which is again a monotonic function of 5 > 0. Notice that by condition (ii), Section [2], we 
have lim^o Y^d) = p. 

Choose e > such that p + e < a/(l + a). Then we have 

F(Ef' N < ^) ^ p(^(Ba/l 2 )>p + e) 
Now we choose <5 such that 

F 0O (5)-p<e/2, 
and let Z be defined by 5 = Bajl\. Then for any I > l we have: 

F{E^' l ' N < B/l 2 ) <: P(V(Ba// 2 ) > p + e) ^ P(V/"(<5) - p > e) 
< P(|F^(5)-F 00 (5)| >e/2) ^ e 



M(5,e)Vi _ 

□ 
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